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of them an exact local Langevin equation, with time-dependent 
coefficients. 
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I. INTRODUCTION 



The master and Langevin equations are the standard approaches to deal 
with irreversible phenomena. However, it is an old problem of physics to obtain 
such macroscopic irreversible equations departing from the reversible micro- 
scopic laws without appealing, a priori, to approximations or extra-dynamical 
hypotheses. 

The simplest systems for which the origin of irreversibility can be studied 
on a microscopic basis are the harmonic linear ones since one is able to 
reduce the Hamiltonian to normal modes, for which the dynamical evolution 
becomes trivial. With this purpose we consider many particle systems with 
Hamiltonians of the form 

N N 

ff = EE (i) 

n=0 m=0 

which conserve the total number of quanta. For this kind of system it can be 
proved by a straightforward calculation that the time evolution of the creation 
operators is given byQ 

bW) = e- im bU0)e im = £ A™(t)&L(0), (2) 

m 

being A nm (t) = (ip m \e~ tht \i() n ^ the transition amplitude between one-particle 
states. 

From Eq. (fj) the exact "Langevin" and "master" equations can be derived. 
This work is devoted to study these equations in order to trace the origin of 



lr rhe meaning of Eq. (|J) is clarified applying to both members of the equality the 
vacuum state. This means that the operator b^t) creates states \ip n (t)) . 
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macroscopic irreversibility. In this route we show, from the reversible evolution 
of the system, which are the simplest, a posteriori, extra-dynamical hypotheses 
needed to understand the irreversible behavior. 



II. THE MASTER EQUATION 

For systems described by a Hamiltonian such as those of Eq. (|I|) we can 
obtain an exact master equation for the mean occupation number (N n ) = 
(b^bn^ , corresponding to the state \ip n ) of the one-particle Hamiltonian h, 
namely 

d(N n (t)) 



J2 W nk(t) (N k (t)} , (3) 



or equivalently in the form of a kinetic balance equation 

d(N n (t)} 



E [Wnm{t) (N m (t)) - W mn (t) (N n {t))} . (4) 

The time-dependent (non-symmetrical) coefficients W nm {t) are given by 

W nk (t) = £ Pnm (t)P m l(t), (5) 

m 

where P nm {t) = |Ai m (£)| 2 is the probability to find a quantum in the state 
\ip m ) at time t if it was in the state \if) n ) at t — 0. 

We can prove Eq. (||) as follows: First, from (fj), it can be checked that| 

(#«(*)> = £ P nm (t) (N m (0)) , (6) 



2 Eq. (^) has a simple interpretation. The mean number of quanta at the level n at 
a given time t can be obtained as a sum of the initial populations of different levels 
times the transition probabilities at t from these levels (including the own n) to the 
level n. 
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provided we assumed a privileged initial condition which represents the absence 
of initial correlations^] 

(bl(O)UO)) = S mn (N m (0)) . (7) 

Finally, as is well known, departing from the solution we can form a differ- 
ential equation by derivation and eliminating the integration constants. That 
is, derivating (J|) with respect to time and inverting this linear system to elim- 
inate (iV m (0)) , we finally obtain (|]). 

III. THE LANGEVIN EQUATION 

Until now the basis {\tp n )} n =o jy is an Y complete set of the one-particle 
space. In this section we consider the case in which the basis {\ip n )} n=0 n 
diagonalizes the unperturbed one-particle Hamiltonian /i (h = h + v), such 
that H = Hq + V can be thought as a set of interacting harmonic oscillators. 
H Q represents the set of uncoupled oscillators and V is a linear interaction 
among them. Splitting the summation in such a way that one oscillator is 
identified as a Brownian particle and the rest as a bosonic reservoir, redefining 
the notation by {\ip ) , l^n) } n =i,.. ., N = {M , kn)} n =i,...,jv > b o = B -> we have 

(h = i) 

N 

H = {IB^B + u n blb n , (8) 

n=l 

and 

N N 

v = (n\v\n) b^b + J2J2 KM^ m > &J>& m + 

n=l m=l 

TV 

£ ((co n \v\n) b\B + (tt\v\uj n ) B^b n ) . (9) 

n=l 



The last hypothesis has the same status as the random phase aproximation. 
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In this case we can derive a generalized form of the Langevin equation for 



the position operator X = ^/ 2 1 Mn [B + B^J of the oscillator with frequency f2 
departing from the exact solution (fj) with n = 0, namely 

X(t) = a(t)X(0) + b(t)^ + f(t), (10) 

1 N 

Ann(t) = a(t) + ib(t), f(t) = -j== £ [A nm {t)b\M + h.c. . 

As in this case we have two constants of integration and a particular solution 
f(t), X(t) satisfies a second-order differential equation like 

x (t) + n 2 (t)x(t) + r(t) x (t) = F{t), (n) 

with an inhomogeneous term given by 

F(t) =f (t) + n 2 (t)f(t) + v(t) f (t), 

which represents the analogue of the stochastic acceleration in the standard 
Langevin equation.f] The unknown coefficients Q 2 (t) and T(t), are the ana- 
logues of the time-dependent square frequency and damping factor, respec- 
tively. They can be easily determined by solving the linear system which 
results in replacing the two independent solutions of the homogeneous equa- 
tions: 

a (t) + n 2 (t)a(t) + T(t) a (t) = 0, 



b(t) + n 2 (t)b(t) + r(t) h(t) = o. 



That is 



q2/£\ __ AgnA nn — A nn A nn ^. . _ Aqq A nn — A nn Ann 
Ann A na -A na Ann Ann A nn -A nn Ann 



However, its deterministic character is obvious from its definition. 



5 



IV. PERTURB ATIVE CALCULATIONS 



The coefficients of the master equation are in general time- dependent. We 
can evaluate the transition probabilities involved in these coefficients using 
time-dependent perturbation theory. The Fermi golden rule gives P nm = S nm + 
T nm t, where 

T nm = 27r \(ip n \ v \ip m )\ 2 5 t (uj n - u m ), for n ^ m, 

r n n = - 2n l(V>n| v \ip m )\ 2 5 t (u n - u m ), 

being 5t(a) = S1 ^" f an approximating of the Dirac delta for very long times. 
Using the last expressions for calculating coefficients W nm (t), neglecting higher 
order terms in the perturbation, we have the time-independent symmetrical 
coefficients W nm (t) = Y.k^nki.bkm - r fcm t) = T nm , in agreement with stan- 
dard results. In the case of the coefficients of the Langevin equation we 
must evaluate up to the second order the survival amplitude of the state \ fl). 
For very long times an exponential contribution dominates its time evolution, 
Ann(t) — e~*( Q+(5n ~*^*, with a frequency shift and a damping factor given by 

sn = (n\v\n) + pJ2 ^ vm \ ^ = 2n\(n\v\n}\ 2 . 

A straightforward calculation shows that in the exponential decay regime the 
coefficients of Eq. (0) are given by Q(t) = Q + 5Q and T(t) — 7 . 

V. CONCLUSIONS AND FURTHER REMARKS 

We have shown that from reversible quantum mechanical laws we can ob- 
tain equations of motion as if they represented stochastic processes. How- 
ever, our equations (0) and QTTJ) describe the exact dynamical evolution of 
the system and in this sense they are equivalent to the Schrodinger and/or 
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Heisenberg equations. We also stress the fact that Eqs. (|3]) and ([0]) are local 
in time in contrast with the standard non-Markovian master and Langevin 
equations. In our case all memory effects are reduced to the knowledge of 
the initial conditions. The time-dependent coefficients of Eqs. (|3|) and (|TT|) 
are uniquely determined for the amplitudes A nm {t) of the one-particle sector. 
Evaluating these amplitudes through perturbation theory, retaining up to the 
second order in the Dyson expansion, we retrieve the standard irreversible 
equations with time-independent coefficients. In Ref. f|| we have considered 
a particular model in which there is not interaction among bath oscillators 
((u> n \ v \u) m ) = 0). This choice allows us to reduce H to normal modes (a set of 
uncoupled harmonic oscillators) , which in the one-particle sector means that h 
is diagonal: h = Y,v=o a v \ a v) (&u\ ■ 111 this case the transition amplitudes can 
be analytically obtained as 

JV 

A nm {t) = WmM M^n) , (13) 

i/=0 

which will allow us to extend the previous analysis beyond the perturbative 
calculations and determine its range of applicability in a further work. 

We thank Leopoldo Garcia-Cohn for fruitful discussions. 
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